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Abstract 

We use a simple model to study the long time fluctuations induced by random pinning on the motion of driven non- 
interacting vortices. We find that vortex motion seen from the co-moving frame is diffusive and anisotropic, with 
velocity dependent diffusion constants. Longitudinal and transverse diffusion constants cross at a characteristic 
velocity where diffusion is isotropic. The diffusion front is elongated in the direction of the drive at low velocities 
and elongated in the transverse direction at large velocities. We find that the mobility in the driven direction is 
always larger than the transverse mobility, and becomes isotropic only in the large velocity limit. 
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1. INTRODUCTION 



2. MODEL 



Understanding the non-equilibrium statistical 
properties of driven particles in disordered media 
is a challenging question relevant to many exper- 
imental situations. A prominent example are the 
moving phases of driven vortex lattices in super- 
conductors [1] . A key feature of these systems is 
that the disorder induces anisotropic response and 
fluctuations which are strongly controlled by the 
velocity [1,2,3]. 

In spite of its relevance to understand situations 
of incoherent or plastic vortex flow the simple case 
of an isolated vortex driven in a d-dimensional ran- 
dom potential with d > 1 has been tackled analiti- 
cally only by perturbation theory [2] , valid at high 
velocities, or by mean field theory [4], valid for d 3> 
1 . In this paper we propose a simple model whose 
long time behaviour can be computed analitically 
at d = 2 for any finite velocity. 



Let us consider the equation of motion, in two 
dimensions, of a driven isolated vortex at zero tem- 
perature, 



V v(t) = F + ^f p ( r - ri ) 



(1) 



where v = dr/dt is the instantaneous velocity of 
the vortex located at r (t) , f is the driving force and 
r\ the friction coefficient. We model the disorder as 
a random arrangement of hard disks with center 
Yi and radius £. Outside the disks the vortex has a 
free motion and inside it feels a pinning force, 



.4 



f p (r) = -^|r9(l-r 2 /e 2 ) 



(2) 



where A p is the amplitude, and O the step func- 
tion. This disorder models a diluted distribution of 
pinning centers separated at a distance d > 2£. In 
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the following we use adimensionalized variables: £ 
is the length unit, A p the energy unit, and A p j(^r\ 
the time unity. 

Above the dcpinning transition the motion of 
the vortex consists in straight segments of free mo- 
tion interrupted by the collisions with the different 
pinning centers. At each collision the vortex is de- 
layed and deflected with respect to the free motion. 
The equation describing the motion of the vortex 
inside the trap centered at r, = is, 



dr 

dt 



(3) 



The collision starts with the vortex at some initial 
position r(0) = r on the border of the trap, r 2 = 
f. The solution of Eq.(3) is 



r(t) = (r - F)e- t + F 



(4) 



After a time interval St the vortex will exit from 
the trap, therefore r 2 (5t) = 1, St > 0. Using this 
condition in equation 4 we obtain the following 
expression for St, 



e~ 5t = 



f°.F-0f p ".FP-(^-f)(f^ 
(f p °) 2 



(5) 



where i p = f p (r ). The displacement Sr induced 
by the collision is then given by, 



Sr = r(St)-r = (r -F)(e' St - 1) 



(6) 



Due to the random distribution of pinning centers, 
the motion can be considered as a random walk 
in the long time limit. The fluctuations of r(t) are 
thus induced by the uncorrelated sequence of colli- 
sions. Assuming identical pinning centers, the ran- 
domness comes exclusively from the random initial 
conditions of each collision, described by Eq. (4) . 

In a long time interval t the vortex collides with 
a large number N c (t) of pinning centers. Since the 
core covers an area ~ Vt£, where V is the mean 
velocity, we get, 



N c (t) w n p £Vt. 



(7) 



By symmetry, the long time displacement is along 
the direction of F, Ar(i) = Ar\\(t)F, and 

Ar,| (t) = Vt « JV c (t)<Jr||> + F(t - N c (t)(St)) (8) 



where (...) denotes an average over a random dis- 
tribution of initial conditions i"o in equations 5 and 
6. From Eq. (8) we get the mean velocity, 



v- 



l-n P ti(S ril )-F(St)] 



(9) 



We can now define longitudinal D\\ and trans- 
verse D± diffusion constants, 

D^{[A n {t)-Vt] 2 )/t (10) 
D ± = (Ari(t)>/t=^<*ri> (11) 

To calculate the longitudinal diffusion constant in 
terms of the single collision displacement we use 
that d — V(FSt + d — Sr»)/F is the random lon- 
gitudinal displacement with respect to the aver- 
age longitudinal motion in a single collision, with 



d = 



_ -1/2 



the longitudinal mean distance between 



pinning centers of two given consecutive collisions. 
We thus get, 



D 



^^([d-jiFSt + d-Sr^)] 2 ) (12) 



We also define the longitudinal mobility as yU|| = 
In terms of single collision quantities we get, 



H\\ = ^{l-Vn P t({5t)+F 



d(St) d(Sr\\) 



dF 



dF 



(13) 



To define the transverse mobility we need to intro- 
duce a small perturbative force f±. The velocity 
induced by this force is, 



v± = Vn P S[(6r ± )-f ± (6t)]+f± 



(14) 



and thus we can define the transverse mobility 
fj,± = [-^-]f ± ->o- In terms of single collision quan- 
tities we get, 



.( d(Sr ± ) 



» ± = l + Vn p £[^ I ^-(5t) 



(15) 



/J.-0 



Finally we can define effective temperatures using 
generalized Einstein relations, 



■D ± /2fx A 
■■D\\/2 n 



(16) 
(17) 



In order to calculate the transport properties de- 
fined above (V, _Dj_, , T e ^ f , and T^ lf ) 
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Fig. 1. (a) VF characteristics of the driven particle. The 
dashcd-line is the free flow solution, (b) Correction to the 
free flow velocity. Lines show the asymptotic forms at small 
and large velocities. 

we need to calculate the first moments of the dis- 
tribution of 5t, d(St)/dF, Sr, 5r\\5t, d(5x±)/df±, 
d{5x\\)/dF by performing simple integrals. 



3. RESULTS 

In Fig. 1(a) we show the Velocity - Force charac- 
teristics (VF) of our model, calculated from Eq.(9). 
The critical dcpinning force is F c = 1. At zero 
temperature, for F < F c , the particle is trapped 
after a transient, and thus V = 0. At low veloc- 
ities, F — > F c + , the VF curve is strongly nonlin- 
ear with V ~ [ra p log((F - F C )/2F C )}' 1 [6] and at 
large velocities free flux flow is approached with 
corrections that scale as F - V ^ V -1 , as shown 
in Fig. 1(b). In Fig. 2(a) we show the longitudi- 
nal and transverse diffusion constants. Dm and D± 
are both non-monotonous functions of the veloc- 
ity V. At small velocity diffusion constants grow 
linearly with V while at large velocity Dm ~ V~ 3 
and D± ~ V~ x . Let us note that Dy and D± cross 
at a characteristic velocity V Q . This crossing means 
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Fig. 2. (a) Diffusion constants in the longitudinal (A) and 
transverse (□) directions. (b) Mobility in the transverse (□) 
and longitudinal (A) directions. Lines show the small and 
large velocity asymptotic forms in both figures. 

that the long time diffusion front changes aspect 
ratio at V ■ For V < V the diffusion front is elon- 
gated in the driven direction, while for V < V is 
elongated in the transverse direction. At V diffu- 
sion is isotropic. Interestingly, the same behavior 
is observed in numerical simulations of interacting 
vortices in two dimensions [3] . 

In Fig. 2(b) we show the longitudinal and trans- 
verse mobilities. We observe that both are velocity 
dependent and approach the free flux response, 
/j, ~ 1, at large velocities. Since the fi\\ is the differ- 
ential resistance, dV/dF, the divergence observed 
at small velocity is a signature of the dcpinning 
transition. near depinning is high since any 
small force can reduce strongly the waiting time 
St ~ 1/V inside the trap. On the contrary, [i± ~ V 
at small velocity. A small transverse force has a 
small effect in the trap time St compared with the 
linear V dependence (Eq.(7)) of the number of 
collisions per unit time. 

In Fig. 3 we show the effective temperatures in 
the longitudinal and transverse directions. Since at 
large velocities the mobilities saturate, the velocity 
dependence of the effective temperatures are dom- 



3 



0.01 




0\ k (a) - 

v- 1 x 

■ 

; \ j 
- 


1Q~ 2 


' — ' ' — 1 — 1 — 111 ^ 

r / \ (b)- 


10" 3 








10" 4 




10~ 5 





1 

V 



Fig. 3. Effective temperatures in the transverse (a) and 
longitudinal (b) directions. Lines indicate the asymptotic 
forms at large velocity. 



inated by the diffusion constants, and thus T e ^ f ~ 
l/V and Tj!j f ~ 1/V 3 . At low velocities the effec- 
tive temperatures are instead strongly determined 
by the velocity dependence of the mobilities. T^ ff 
reaches a maximum at a intermediate velocity and 
decreases quickly to zero with decreasing V. On 
the other hand, T^ f saturates at a finite value as 
the velocity vanishes, due to the linear dependence 
with velocity of both the diffusion constant and the 
mobility. 

It is worth noting that only at large veloci- 
ties and in the transverse direction, the effective 
temperature is found to be identical, except from 
numerical factors, to the shaking temperature 
[2]. This confirms that the shaking temperature 
is equivalent to an effective temperature defined 
from a generalized fluctuation relation only in the 
limit of large velocity of non-interacting vortices 
[3] . It is important to point out here that the ther- 
modynamic nature of the effective temperatures 
of this model is still unclear, since the system is 
non-interacting and strongly driven [5] . 



4. CONCLUSIONS 



We have studied the pinning-induced anisotropic 
diffusion of driven non-interacting vortices. We 
find that the diffusion front is elongated in the 
direction of the driving force at low velocities 
and in the transverse direction at large velocities. 
This implies the existence of isotropic diffusion 
at a characteristic velocity V . The analysis of 
the anisotropic low frequency voltage noise in su- 
perconductors could be a possible experimental 
probe of this result, since diffusion constants and 
velocity fluctuations are related by generalized 
Green-Kubo relations [3] . 

Even if the depinning transition does depend 
on the peculiarities of our model, we find that 
the main features of the long time fluctuations we 
obtain at intermediate and large velocities are in 
agreement with perturbation theory predictions [2] 
and with numerical simulations of non interacting 
vortices for the non-simplified model [6] . Further- 
more, it is easy to show that our model can be 
solved at any dimension d > 1 and generalized to 
more complicated short-range random potentials 
[6]. 
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